Simulation Studies on the Stability 
of the Vortex-Glass Order 
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The stability of the three-dimensional vortex-glass order in random type-II superconductors with 
point disorder is investigated by equilibrium Monte Carlo simulations based on a lattice X Y model 
with a uniform field threading the system. It is found that the vortex-glass order, which stably 
exists in the absence of screening, is destroyed by the screenng effect, corroborating the previous 
finding based on the spatially isotropic gauge-glass model. Estimated critical exponents, however, 
deviate considerably from the values reported for the gauge-glass model. 



Due to the enhanced effect of fluctuations, the prob- 
lem of the phase diagram of high-T c superconductors in 
applied magnetic fields is highly nontrivial and has at- 
tracted much interest recently. For random type-II su- 
perconductors with point disorder, possible existence of 
an equilibrium thermodynamic phase called the vortex- 
glass (VG) phase, where the vortex was pinned on a long 
length scale by randomly distributed point-pinning cen- 
ters, was proposed In such a VG state, the phase of 
the condensate wavefunction is frozen in time but ran- 
domly in space, with a vanishing linear resistivity pi,. 
It is a truly superconducting state separated from the 
vortex-liquid phase with a nonzero pl via a continuous 
VG transition. 

This proposal was supported by subsequent experi- 
ments. In particular, transport measurements on films 
H and twinned single crystals (|] gave evidence for the 
occurrence of a continuous transition into the glassy su- 
perconducting state. It should be noticed, however, that 
these samples often contain extended defects with corre- 
lated disorder, such as grain boundaries, twins and dis- 
locations. Since these extended defects could pinn the 
vortex more efficiently than point defects, the possibility 
still remains that extended defects play a crucial role in 
the experimentally observed "VG transitions", and the 
sample only with point defects behaves differently [4-7]. 

Stability of the hypothetical VG state was also studied 
by numerical simualtions. Here, it is essential to obtain 
the data of appropriate thermodynamic quantities in true 
equilibrium and to carefully analyze its size dependence. 
So far, such calculations have been limited almost exclu- 
sively to a highly simplified model called the gauge-glass 
model. Previous simulations on the three-dimensional 
(3D) gauge-glass model have indicated that, while the 
stable VG phase exists in the absence of screening ||^| , 
the finite screening effect inherent to real superconduc- 
tors eventually destabiizes it pQ,Q|. 



Meanwhile, it has been recognized that the gauge-glass 
model has some obvious drawbacks ||. First, it is a spa- 
tially isotropic model without a net field threading the 
system, in contrast to the reality. Second, the source of 
quenched randomness is artificial. The gauge-glass model 
is a random flux model where the qeunched randomness 
appears in the phase factor assoicated with the flux. In 



reality, the quenched component of the flux is uniform, 
nothing but the external field, and the quenched ran- 
domness occurs in the superconducting coupling. It re- 
mains unclear whether these simplifications underlying 
the gauge-glass model really unaffect the basic physics of 
the VG ordering in 3D. 

The purpose of the presnet letter is to introduce a 
model in which the above limitations of the gauge-glass 
model are cured, and to examine by extensibe Monte 
Carlo (MC) simulations the nature of the 3D VG order- 
ing with and without screening beyond the gauge-glass 
model. It is found that, as in the gauge-glass model, the 
VG phase is stable in the absence of screening but the 
finite screening effect destabilizes it. 

We consider the dimensionless Hamiltonian[12,13], 
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where J is the typical coupling stregth, 9i is the phase 
of the condensate at the i-th site of a simple cubic lat- 
tice, A is the fluctuating gauge potential at each link of 
the lattice, the lattice curl V x A is the directed sum of 
Ay's around a plaquette with Aji — —A^, and Ao is the 
bare penetration depth in units of lattice spacing. <fr ex t 
is an exernal flux threading the elementary plaqutte p, 
which is equal to h if the plaquette is on the a:y-plane and 
zero otherwise, i.e., a uniform field is applied along the 
z-direction. The first sum in (1) is taken over all nearest- 
neighbor pairs, while the second sum over all elementary 
plaquettes. Fluctuating variables to be summed over are 
the phase variables, Oi, at each site and the gauge vari- 
ables, Aij, at each link. In order to allow for the flux 
penetration into the system, we impose free boundary 
conditions in all directions for both 9i and Aij [12,13]. 
Quenched rondomness occurs only in the superconduct- 
ing coupling Jij which is assumed to be an independent 
random variable uniformly distributed between [0,2]. We 
stress that the aforementioned drawbacks of the gauge- 
glass model have been cured now: The present model has 
a uniform field threading the system and the quenched 
randomness occurs in the superconducting coupling, not 
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in the gauge field. 

In addition to the global U(l) gauge symmetry, the 
Hamiltonian (1) has a local gauge symemtry, i.e., the in- 
variance under the local transformation 9i — > 6i + A and 
Ai^i+s — * A-i.i+s + A for an arbirtary site i (Ai^+g's are 
all link variables emanating from the site i). We adopt 
the Coulomb gauge, imposing the condition, divA = 
Ai t i + s = 0, at every site i. In the limit of vanish- 
ing screening Ao — > oo, the link variable Aij is quenched 
to the external-field value, and the fluctuating variable 
becomes the phase variable 9i only. 

Simulation is performed based on the exchange MC 
method where the systems at neighboring temperatures 
are occasionally exchanged |14|. We run two indepen- 
dent sequences of systems (replica 1 and 2) in paral- 
lel, and compute a complex overlap q between the lo- 
cal superconducting order parameters of the two replicas 
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where the summation is taken over all N = L 3 sites. In 
terms of the overlap q, the Binder ratio is calculated by 
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where (• ■ •} represents the thermal average and [• ■ •] rep- 
resents the average over bond disorder. Note that, thanks 
to the Coulomb-gauge condition, the superconducting or- 
der parameter, which is originally not local-gauge invari- 
ant, becomes a nontrivial quantity. 

We deal mainly with two cases; [I] no screening cor- 
resonding to Ao = oo, and [II] finite screening corre- 
sponding to Ao = 2, whereas some data are taken for 
the case of stronger screening corresponding to Ao = 1. 
In either case, we fix the field intensity to h = 1 which 
corresponds to / = l/(2n) ~ 0.159 flux quanta per pla- 
quette. We have chosen the fractional value of / to 
avoid the commensurability effect associated with the 
vortex-lattice formation. The lattice sizes studied are 
L = 6,8,10,12 and 16 (A = oo), and L = 6,8,10,12 
(Ao = 2). Equilibration is checked by monitoring the 
stability of the results against at least three-times longer 
runs for a subset of samples. Sample average is taken 
over 300 [L = 6), 200-300 (L = 8), 120 (L = 10), 75-150 
(L = 12) and 136 (L = 16) independent bond realiza- 
tions. 

We begin with the case of no screening (Ao = oo). 
The size and temperature dependence of the calcualted 
Binder ratio is shown in Fig. 1(a). As can be seen from 
Fig. 1(a), g(L) for different L > 8 tends to merge, or 
weakly cross, at T/J = 0.68 ± 0.02, indicating that the 
VG transition occurs at a finite temperature in the ab- 
sence of screening. Observed near marginal behavior sug- 
gests that D = 3 is close to the lower critical dimension. 




FIG. 1. Temperature and size dependence of (a) the bulk 
Binder ratio, Eq.(3), and of (b) the transverse Binder ratio, 
Eq.(5), of the model without screening (Ao = oo). The insets 
represent magnified views of the transition region. 

In the present model, as well as in reality, the nature of 
fluctuations along the field (longitudinal direction) and 
perpendicular to the field (transverse direction) could dif- 
fer. An extreme possibility here may be that the VG or- 
der occurs only in some spatial component, say, in the 
transverse component, keeping the other (longitudinal) 
component disordered Indeed, the possibility of such 
"two-dimensional" or purely "transverse" vortex order 
has been discussed in the literature as a "decoupling" 
transition |l5[ |. In order to probe such exotic possibil- 
ity, we introduce a transverse Binder ratio in terms of 
the layer-overlap q' k defined for the fc-th xy-layer of the 
lattice by, 
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When the VG order occurs in each layer, gtrans(L — > oo) 
should be nonzero. In particular, if the purely trans- 
verse VG order is to occur as a consequence of the 
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layer-decoupling, gtrans(L — > oo) should stay finite while 
g[L — > oo) vanishes. 

The calculated gtrans(L) is shown in Fig. 1(b). As can 
be seen from Fig. 1(b), gtrans(L) exhibits behavior quite 
similar to <?(£), revealing a merging or a weak crossing 
at T/J = 0.67 ± 0.02. This indicates that the present 
model exhibits only a single bulk VG transition where 
both the transverse and longitudinal components order 
simulataneously. 




FIG. 2. Temperature and size dependence of (a) the bulk 
Binder ratio and of (b) the transverse Binder ratio, of the 
model with screening (Ao = 2). The insets represent magni- 
fied views. 

In the case of hnite screening Ao = 2, as shown in 
Figs. 2(a) and 2(b), g(L) and gtrans{L) constantly de- 
crease with increasing L at all temperatures studied, sug- 
gesting that a finite-temperature VG transition is absent 
in the presence of screening. Similar behavior is observed 
also for the case of stronger screening, Ao = 1, and it can 
be concluded that the screening effect destabilizes the 
VG order at finite temperature. Thus, concerning the 
presence or absence of the VG order, the present model 
yields the same answer as the spatially isotropic gauge- 
glass model. 

Next, we turn to the critical properties of the model. 
For a finite-temperature VG transition for Aq = oo, we 



estimate the correlation-legth exponent v = 2.2 ± 0.4 via 
the finite-size scaling analysis of g(L) (see Fig. 3). Then, 
from the order parameter = [< q 2 >], the critical- 
point-decay exponent is determined to be 77 = —0.5 ±0.2. 
Likewise, from the decay of the autocorrelation function 
of the superconducting order parameter at T = T g [ [l6|| , 
the dynamical exponent z is estimated to be z = 3.3±0.5. 
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FIG. 3. Finite-size scaling plot of the bulk Binder ratio of 
the model without screening (Ao = 00). 

In order to examine the possibility of anisotropic scal- 
ing where the longitudinal and transverse correlation 
lengths have different v values, we perform the scal- 
ing analysis of gtrans and q trans as compared with g 
and q^ |l^]: If the scaling is really anisotropic, there 
is no reason to expect that the exponents determined 
from these transverse quantities coincide with those de- 
termined from the bulk quantities. We get v ~ 1.8 and 

(2) 

77 ~ -0.5 from 

Qtrans 

and qt rans , which agree within the 
errors with v and 77 determined above from g and q( 2 \ 
We found no evidence of anisotropic scaling, although the 
occurrence of small anisotropy cannot be ruled out from 
the present data. 

Concerning the T g = transition of Ao = 2, we get 

v ~ 4.0 and v ~ 3.5 from g and q( 2 \ and v ~ 4.0 and 

(2) 

v ~ 3.0 from gtrans and ql rans , which again agree within 
the errors. (Note that r\ should be equal to -1 for a T = 
transition with nondegenerate ground state as expected 
in the present model.) Hence, the scaling also appears to 
be isotropic with v = 3.5 ± 1.0. 

The obtained exponents are summarized in Table 1 and 
are compared with the values reported for the spatially 
isotropic gauge-glass model with and without screening. 
Apparently, there exists a significant deviation between 
the two results. In particular, v of the present model 
appears to be significantly larger than v of the gauge- 
glass model which might suggest that the present 
model lies in a universality class different from that of 
the gauge-glass model. 

Finally, we wish to discuss the experimental implica- 
tion of our results. The present study suggests, in accord 
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TABLE 1 Critical exponents of the present model com- 
pared with the values reported for the gauge-glass model 
for both cases of Ao = oo (no screening) and Ao < oo 
(finite screening). 





model 


v rj z 


Ao = oo 
(T g > 0) 


present 


2.2(4) -0.5(2) 3.3(5) 


gauge glass 




Ao < oo 
(T g = 0) 


present 


3.5(10)_ -1 


gauge glass 


1.05(10) 1 10 1 
1.05(3) [11] 



with the previous studies based on the gauge-glass model, 
that in random type-II superconductors with point dis- 
order there should be no VG phase at finite temperature 
in the strict sense. As discussed, experimental data on 
films and twinned crystals supporting the occurrence of a 
thermodynamic transition into the truly superconducting 
glassy state might well reflect the properties associated 
with extended detects. In this connection, the proper- 
ties of a sample exclusively containing point defects is 
of great interest [4-6] . Recently, Petrean et al measured 
the transport properties of such sample, untwinned, pro- 
ton irradiated YBCO JXSj ] . These authors observed an 
Ohmic behavior at all temperature studied, but found 
that the linear resistivity appeared to vanish at a finite T* 
as pl ~ (T— T*Y with a universal exponent s ~ 5.3±0.7. 
Since the non-Ohmic region could not directly be reached 
in these measurements, it is not clear at the present stage 
whether the non-Ohmic behavior really sets in below the 
apparent T* deduced from the high-temperature Ohmic 
regime. An another possibility, which is entirely consis- 
tent with the present result, may be that the power-law 
decrease of Pl eventually breaks down at some temper- 
ature close to T* , yielding a small but nonzero pz, even 
at T < T*. If this is the case, the universal critical be- 
havior oi pl observed above T* should be governed by 
the Ao = oo VG fixed point [1,20] which, however, should 
eventually be unstable against the screening effect. In- 
deed, our present estimate of s = (z — l)v 5.1 for the 
Ao = oo transition is close to the experimental value of 
Ref . |h| . It might be interesting to experimentally deter- 
mine the exponents v, z and r\ separately, as well as to go 
further down to lower temperatures in order to examine 
whether the non-Ohmic behavior really sets in below T* . 

In summary, we have introduced a VG model which 
possesses a uniform field and cures the limitations of the 
gauge-glass model. Eextensive simulations show that, 
while the stable vortex-glass phase occurs in the absence 
of screening, it is eventually destroyed by the screening 
effect. Critical exponents associated the VG transitions 
appear to differ from those reported for the gauge-glass 
model, suggesting that real VG transitions may lie in a 
different universality class. 



The numerical calculation has been performed on the 
FACOM VPP500 at the supercomputer center, ISSP, 
University of Tokyo. The author is thankful to R. Ikeda 
and S. Okuma for useful discussion. 
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